Aspects of the confinement mechanism in Coulomb-gauge QCD 
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Phenomenological consequences of the infrared singular, instantaneous part of the gluon propa- 
gator in Coulomb gauge are investigated. The corresponding quark Dyson-Schwinger equation is 
solved, neglecting retardation and transverse gluons and regulating the resulting infrared singular- 
ities. While the quark propagator vanishes as the infrared regulator goes to zero, the frequency 
integral over the quark propagator stays finite and well-defined. Solutions of the homogeneous 
Bethe-Salpeter equation for the pseudoscalar and vector mesons as well as for scalar and axial- 
vector diquarks are obtained. In the limit of a vanishing infrared regulator the diquark masses 
diverge, while meson properties and diquark radii remain finite and well-defined. These features are 
interpreted with respect to the resulting aspects of confinement for colored quark-quark correlations. 
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The substructure of the nucleon has been determined 
to an enormous precision leaving no doubt that the par- 
ton picture emerges from quarks and gluons, the ele- 
mentary fields of Quantum Chromodynamics (QCD). Al- 
though they are the "elementary particles" of strong in- 
teractions, quarks and gluons have never been detected 
outside hadrons. This phenomenon is called confinement. 
Despite its importance for particle physics and for an 
axiomatic approach to quantum field theory our under- 
standing of confinement is far from being satisfactory. 

In this letter we concentrate on certain aspects of con- 
finement for colored composite states. We start from the 
commonly accepted Wilson criterion and an inequal- 
ity between the gauge-invariant quark-antiquark poten- 
tial Vw(R) and the color-Coulomb potential Vc(x) 0. 
The latter quantity is the instantaneous part of the time- 
time component of the gluon propagator in Coulomb 
gauge: Doo(x,t) oc Vc(x)S(t)+ non-inst. terms. In 
Ref. it was shown that if Vw{R) is confining, i.e. if 
limii-xx, Vw(R) —> 00 ; then also |Vc7(x)| is confining. 
This was confirmed in an SU(2) lattice calculation |3| 
where it was found that — Vc(x) rises linearly with R = 
\x\. However, the corresponding string tension, cr c , was 
extracted to be several times the asymptotic one 1 . 

A well-suited formalism for the study of composite or 
bound states of quarks is the Dyson-Schwinger/Bethe- 
Salpeter approach [f|. While corresponding investiga- 
tions in Coulomb gauge, e.g. 0,0, predate those based 
on model studies employing Landau-gauge QCD Green 
functions, the latter have been much more numerous 
and the corresponding studies explore a large number 



1 If the same holds for the physical case of three colors one infers 
u c as 600 . . . 750 MeV from the generally used value u c rj 440 
MeV. Note, however, that this increase is not sufficient to resolve 
the problem of a too small value of the pion decay constant Q , 
when only a confining potential is used and non-instantaneous 
interactions, in particular transverse gluons, are neglected. 



of hadron observables, see e. g. Refs. [&i @, Hfl Ull IL2I 

and references therein. Note that in Landau-gauge QCD 
the structure of the quark-gluon vertex 0] is an issue 
of current debate due to its importance for the quark 
propagator [lij . 

In this letter we report on a study of mesons and 
two-quark composite states employing the color-Coulomb 
potential Vc(x) and thus some of the basic features of 
Coulomb-gauge QCD. We build on investigations of the 
gluon propagator 1 fill and the dynamical breaking of chi- 
ral symmetry ,ll7j in Green-functio n ap proaches and 
related results of lattice calculations 0, Il8j| . Our focus is 
the realization of confinement for quarks and two-quark 
composite states ("diquarks"). 

First, we briefly review the quark Dyson-Schwinger 
(gap) and bound-state Bethe-Salpeter equations. All cal- 
culations are performed in Minkowski space. The QCD 
gap equation determines the quark self-energy due to glu- 
ons. It is of the form 



(1) 



where S(p) is the renormalized dressed quark propaga- 
tor, m the current-quark mass, and £(p) is the quark self 
energy. A quark-antiquark bound state is described by 
the Bethe-Salpeter equation (BSE), which in its homoge- 
neous form is written as (for simplicity we neglect Dirac, 
flavor, and color indices) 



re?,?) 



d 4 k K(q, k, P) S(k+) T(P, k) S{k_) , (2) 



where P and q are the quark-antiquark pair's to- 
tal and relative four-momenta, T(P, q) is the bound 
state's Bethe-Salpeter amplitude (BSA), k± = k ± P/2 
are the individual quark- and antiquark-momenta, and 
K(q, k, P) is the quark-antiquark scattering kernel. Note 
that the result of Eq. (JTJ appears as input in Eq. . 

The quark self-energy in Eq. (JJJ is a functional of the 
quark and gluon propagators and the quark-gluon vertex; 
a self-consistent solution would require to simultaneously 
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FIG. 1: The pion and rho masses as functions of the current- 
quark mass. All quantities are given in appropriate units of 
y/ol, <J C being the Coulomb string tension (see text). 



FIG. 2: The quark mass function M(q 2 ) for four values of the 
infrared regulator fim. All quantities are given in appropriate 
units of -v/av. 



solve the Dyson-Schwinger equations for these functions 
and the quark-antiquark scattering kernel from Eq. (J2J 
as well. However, these equations again involve higher 
Green functions and therefore a truncation of this infinite 
coupled system of integral equations is necessary. 

In the present study we use Coulomb gauge together 
with an instantaneous approximation and neglect the ef- 
fects of transverse gluons. These approximations simplify 
the technical challenges involved in concrete calculations. 
On the other hand, some component of the physics con- 
tained in the system is lost. The results are qualita- 
tively, but not quantitatively significant. We therefore re- 
frain from using physical dimensions, but instead present 
the quantities in all graphs in appropriate units of the 
Coulomb string tension a c . The reason for the qualita- 
tive reliability of the calculations is that the underlying 
symmetries of the theory are incorporated in the model 
via Slavnov- Taylor or Ward-Takahashi identities. One 
important example is the axial-vector Ward-Takahashi 
identity, which is used to ensure that the kernels of the 
gap and Bethe-Salpeter equations for pseudoscalar states 
are related in such a way that chiral symmetry and its 
dynamical breaking are respected by the truncation. In 
particular this leads to the correct behavior of the pion 
mass as a function of the current-quark mass in the chi- 
ral limit. This behavior is shown in Fig. ^ (Note: In 
the following we will in general present results for the 
chiral limit. The results for finite current-quark mass are 
analogous.) In this way one can reliably make qualitative 
statements about hadrons and their properties; however, 
it is still important to investigate the contributions from 
retardation effects and transverse gluons, and such efforts 
are currently made. 

In our model the quark self energy S(p) in Eq. (JTJ takes 
the form 

E(p) =(7/6* J j^Vc (k) 70 S{q) 7o , (3) 
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FIG. 3: The masses of the p as well as the scalar (SD) and 
axial-vector (AD) diquarks as functions of the infrared regu- 
lator fim. The mass of the it is identically zero for all values 
of fj,XR an d therefore not shown in the graph. All quantities 
are given in appropriate units of y/aZ. 

where C f = (N 2 - l)/(2iV c ) = 4/3 and k = p-q. In 
the following we will use p to denote p = The go- 
integration in Eq. J2J| can be performed easily. One makes 
the Ansatz S~ 1 (p) := — «(7oPo — 7 ■ P C(p) — B(p)) and 
obtains two coupled integral equations for the functions 
B{p) and C{p) 

B<p)= m+1 ^J* q Vo(k)%$ (4) 

C(p) = 1 + ( d\ V c {k)p-q , (5) 

where p = p/p, m is the current-quark mass, uj(p) := 
v / M 2 (p)+p 2 , and M{q) := B{q)/C{q) is the quark 
"mass function" . Its infrared behavior is a result of dy- 
namical chiral symmetry breaking and can be used to de- 
fine a constituent-quark mass; we have plotted the mass 
function M as a function of q 2 in Fig. (details of this 
figure will be specified below). 

The same approximations and conventions are used in 
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FIG. 4: Pion Bethe-Salpeter amplitude components g and h 
as functions of the infrared regulator [iir,. For convenience, 
the amplitudes are normalized such that h(0) — 1. All quan- 
tities are given in appropriate units of ^foZ. 



the BSE. For pseudoscalar mesons in our model (and 
correspondingly scalar diquarks) the BSA can be char- 
acterized in terms of two scalar functions h(j>) and g(p), 
which essentially are the coefficients of the pseudoscalar 
and axial-vector structures in the BSA. For details, see 
Ref. [H. The BSE, Eq. ©, in terms of h(p) and g(p) in 
our model becomes 



h{p)u{p) = JL Jd 3 qV c (k) 

o 

1TJ. 

g{p) [u(p] 



h(q) 



(G) 



^fd 3 qV c (k) 



h(p) + 

M(p)M(q)+p-q 
Q(p)u)(q) 



9(1), (7) 



where is the bound state's (e. g. the pion's) yet un- 
known mass and uj{p) = C(p) oj{p). 

For vector mesons (and correspondingly axial-vector 
diquarks) the BSA has four linearly independent am- 
plitudes. The construction of the four coupled integral 
equations corresponding to the BSE is analogous to the 
pseudoscalar case. 

The Coulomb-gluon part Vc of the interaction in 
Eqs. (©, © and © is chosen to be 



V c (k) = 



(fc 2 ) 2 



(8) 



where a c is the Coulomb string tension. Obviously, Vc(k) 
is infrared singular. It is regulated by a parameter piR 
such that the momentum dependence is modified to 



V c (k) = 



(fc 2 ) 2 (fc 2 +Mi 2 R ) 2 



(9) 



In this fashion all quantities and observables become 
/j,iR-dependent and one obtains the final result for some 
f(pm) by taking the limit / = lim MlR _ +0 /(Mm)- This is 
illustrated for the quark mass function in Fig. El M(p 2 ) 



FIG. 5: Same as Fig. [I]for the scalar diquark. 



is plotted for different values of /Uir and it is clear that 
the curves converge onto a final result for pm — > 0. 

In order to check the UV behavior one can use a 
Richardson potential |l9j . which has the momentum de- 
pendence V c {k) ~ l/(fc 2 ln(l + fc 2 /A 2 )). The advan- 
tage of our choice for Vc is that the angular integra- 
tion required to solve Eqs. (0J, JSJ can be performed 
analytically. We have checked that the qualitative re- 
sults presented in this paper can be reproduced with the 
Richardson potential. Details of this approach and its 
UV renormalization will be published elsewhere. 

The homogeneous BSE in Eq. J2| is solved by introduc- 
ing an eigenvalue X(P 2 = M 2 ) with M the bound-state 
mass. One then finds M such that A = 1 (for mesons) 
and A = 2 (for diquarks). For details, see e. g. |l2j. 

The curve for A(P 2 ) gets less inclined with smaller 
values of the infrared regulator pm, 1 and its intersection 
point with A = 1 stabilizes in the limit /xir. — > 0. As a 
consequence, while the meson mass is stable, the mass 
eigenvalue for the corresponding diquark state (corre- 
sponding to A(M) = 2) increases like 1//!ir, ultimately 
completely removing these states from the physical spec- 
trum. We have illustrated these effects in Fig.|3|for values 
of 1CT 4 < p m < lCT 2 . 

Note: for N c — 2 diquarks correspond to baryons. In 
particular, in ladder approximation the respective color 
factors for meson and baryon BSEs are identical. There- 
fore the properties of the scalar (axialvector) baryon are 
identical to those of the pion (g meson). For N c > 3 the 
ratio of quark-quark to quark-antiquark color factors in- 
creases like N c — 1; this means that the argument given 
above is also valid in the large- N c limit. 

We studied the BSAs as /im — > 0: the results for g 
and h I© are presented in Figs. 0] and |3] for the pion 
and scalar diquark, respectively. For convenience, the 
normalization of the amplitudes has been chosen such 
that h(0) = 1. We note, however, that IR-cancellations 
appearing in the pion case lead to a stable h as well as 
ratio of g/h, which is not the case (as one would naively 
expect) in the diquark case: there g/h ~ /iiR — > and 
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FIG. 6: Charge radii for tt meson as well as scalar diquark 
(SD) as functions of the infrared regulator nm- All quantities 
are given in appropriate units of ^/a^. 

h ~ W^iR - Still, one can investigate the charge radii 
for both meson and diquark states by requesting that 
the electromagnetic form factor at the origin yields the 
bound-state charge, which gives finite results in the limit 
fJ-m. — > 0. Plots of the pion and scalar diquark charge 
radii are shown in Fig. The results for vector-meson 
and axial-vector-diquark amplitudes are analogous. 

We have performed a study of pseudoscalar- and 
vector-meson states and their corresponding diquark 
partners in a simple model of Coulomb-gauge QCD in the 
context of Dyson-Schwinger equations, which allows for 
obtaining reliable qualitative information about hadrons. 
The infrared singularities in the integrands are regulated 
by the scale /im such that final results are obtained in 
the limit ^ir — > 0. In this limit the masses and charge 
radii for the mesons are stable; for their diquark partners 
only the masses diverge like 1//zir, while the charge radii 
do not. Thus the diquarks are removed from the physical 
spectrum reflecting confinement of colored quark-quark 
correlations. Nevertheless they possess a well-defined 
size. This adds to the motivation of nucleon studies in a 
covariant quark-diquark picture [Tot Hl| . 
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